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Abstract. We prove compactness and hence existence for solutions to a 
class of non linear transport equations. The corresponding models com- 
bine the features of linear transport equations and scalar conservation laws. 
We introduce a new method which gives quantitative compactness estimates 
compatible with both frameworks. 

1 Introduction 

Recent developments in the modeling of various complex transport phenom- 
ena (from bacteria to pedestrians' flows) have produced new and challenging 
equations. In particular those models have a very different behaviour from 
the usual fluid dynamics when the density is locally high, usually as a con- 
sequence of a strict bound on the maximum number of indivivuals that one 
can have at a given point. 

The mathematical theory for well posedness and particularly existence is 
still however lacunary for those equations. The aim of this article is thus to 
provide a unified framework for a general class of conservation laws, including 
many of these recent additions. More precisely, we study equations of the 
form, 

d t n( l t,x)+dxv(a(t,x)f(n(t,x))) = 0, t G K + , x G M d , (1.1) 

where n usually represents a density of individuals. / G W 1,00 (EL, R) is a 
given function which takes local, non linear effects into account. A typical 
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example for / is the logistic f(n) = nil — n/n) + , which limits the velocity 
of individuals when their density is too high thus ensuring that the density 
never exceeds a critical value n. The field provides the 

direction for the movement of individuals. 

Depending on the exact model, a can either be given or be related to 
n. Many such models have been introduced in the past few years in various 
contexts from chemotaxis for cells and bacteria to pedestrian flow models. 
We only give here a few such examples. 

Typically a incorporates some non local effects on the density such as 
with a convolution a = K * n or a Poisson eq. 

a(t,x) = -V x <j)(t,x), -A x (f)(t,x) = g(n(t,x)), (1.2) 

where g is another given function of n. Such a model was introduced in two 
dimensions and in the context of swarming in |33j . The same kind of models 
was studied in [B] and [T5] for chemotaxis (and typically for g{n) = n). 

More complicated relations between a and n are possible, for instance a 
Hamilton- Jacob i equation as in [17] 

a(t,x) = -V x (f>(t,x), -A x .0(t,x) + a|V0| 2 = g(n(t,x)), (1.3) 

with a > (possibly vanishing) and again g a given non linear function. 

Eq. fll.ip can be seen as a hybrid model, combining features of usual 
linear transport equation and scalar conservation laws. 

Let us briefly discuss the main difficulty in obtaining existence of distri- 
butional solutions to fll.ip . With reasonable assumptions (like / ~ n(l — n)), 
it is easy to show that the density n is bounded in every L p spaces. However 
contrary to linear transport equations, a bound on n is not enough to pass 
to the limit in the nonlinear term f(n) (or g(n) if (11.21) or (jl.3p is used). 

With Eq. (TO]) or flEJ} and n G L 1 n L°°, one can easily get a G IV 1 * for 
any 1 < p < oo. From that one may obtain compactness on a in L} oc . 

Hence as a non linear model, the main difficulty in obtaining existence of 
solutions to fll.H) is to prove compactness for the density n. Below we briefly 
indicate why the usual methods for conservation laws do not work in this 
setting (see [13] or [32] for more on conservation laws). 

When a is regular enough (Lipschitz more precisely), then the usual 
method of compactness for scalar conservation laws work and one can for 
example show propagation of BV bounds on n. Unfortunately this Lipschitz 



2 



bound does not hold here in general (only W 1,p , p < oo as explained above). 
Such BV bounds on n can in fact only be propagated for short times (see [B] 
for instance). 

For scalar conservation laws, another way to obtain compactness is either 
by compensated compactness or other regularizing effects. However in di- 
mension larger than 1, those cannot be used as the flux cannot be genuinely 
non linear (it is in only one direction, the one given by a). The 1-dimensional 
case is quite particular (not only in this respect) and many well posedness 
results have already been obtained (see for instance [TT]). 

As far as we know, [15] is for the moment the only result showing ex- 
istence to an equation like (II. ip over any time interval and any dimension. 
The authors use a kinetic formulation of (II. ip . which simply generalizes the 
kinetic formulation of scalar conservation laws introduced in [25] (see also 
[26] and [31]). A rigidity property inherent to the kinetic formulation then 
provides compactness. However a precise connection between a and n is 
needed; more precisely the result is obtained only for the case of (ll.2j) (with 
g = Id though it can obviously be extended to any g suitably regular). 

We conclude this brief summary of the various techniques already in use 
by mentioning gradient flows. In the context of the non linear model ( II. ip . 
the theory is essentially still in development. It requires a lot of structure 
on the equations and that essentially means for the moment Eq. (II. 2p with 
g = Id (any generalization to non linear g would be problematic). We refer 
in particular to [IS] where the right metric for the problem and its properties 
are introduced and studied. 

Gradient flows techniques were also used in [28J for a related problem. 
In that case the corresponding transport is linear but associated with a 
constraint on the maximal density. In the framework of (II. ip that would 
correspond to /(£) = £ I^<i. 

Let us now formulate the main results of the paper. Consider a vanishing 
viscosity approximation 

d t n £ (t, x) + div (a E (t, x) f(n £ (t, x))) - e 2 A x n e = 0, teR+, xE R d , ^ 
n £ {t = 0, x) = n®(x). 

Instead of assuming a precise form or relation between a £ and n e , we make 
very general assumptions on a e . Assume that on [0, T] 

3p>l, sup sup \\a E (t, ■)\\wi>p(K d ) < °°, (1.5) 

e i€[0, T] 
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SUp ||div a! O e || I , o([o ) T]xK") < oo. 



(1.6) 



£ 



As for linear transport equation, an additional condition is needed on the 
divergence to obtain compactness. In order to be compatible with (11.21) or 
(|1.3p . we assume 



Then one can prove 

Theorem 1.1 Assume (jl.5p . (11.61) . (11.71) . that a £ is compact in LP, that n° £ 
is uniformly bounded in L 1 fl L°°(IR d ) and is compact in L 1 (R d ). Then the 
solution n £ (t,x) to (II. 4p is compact in Lj oc ([0, T] x R d ). 

This in particular implies existence results like 

Corollary 1.1 Assume that f E W 1 ' 00 , g E C 2 , /(0) = g(0) = and that 
/(£)#(£) > — C |Cl for some given constant C. Letn E L 1 nL 00 (R d ), a > 0, 
t/ien 3n G L°°([0, T], L 1 fl L°°(M. d )) solution in the sense of distribution 
to (II. ip with (II. 2p . Moreover n is an entropy solution to (11.11) in the usual 
sense that V0 E C 2 convex, 3q E C 1 s.t. 

d t (4>(n(t, x))) + div x (a{t, x) q(n(t, x))) + (0'(n) f(n) - q(n)) div x a < 0. 



Note that this is just one example of possible results, it can for instance 
easily be generalized to (II. 3p under corresponding assumptions. Once a is 
given and in W 1,p the uniqueness of the entropy solution to (11.11) is actually 
not very difficult. However uniqueness for a coupled system like (ll.ll) - (ll.3p 
is more delicate and left open here. 

To prove Th. 11.11 we develop a new method which is a sort of quantified 
version of the theory of renormalized solution and compatible with the usual 
L l contractivity argument for scalar laws. 



Renormalized solutions were introduced in [18] to prove uniqueness to 

solutions of linear transport equations 




r e (x) -r e (y)\ < C\n £ (t,x) - n £ (t,y)\. 



d £ + r £ with d £ compact and 



(1.7) 




d t n + div (a n) = 0. 
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The compactness of a sequence of bounded solutions is obtained as a con- 
sequence of the uniqueness (by proving for instance that w — lim fc n^ = 
(w — linn„. n*;) 2 ). The theory was developped in [18] for a G W 1 ' 1 with 
diva G L°°. It was later extended to a G BV, first for the particular case 
of kinetic equations in jl] (see also [8] for the kinetic case with less than one 
derivative on a). The general case was dealt with in |TJ (see also [EE]). For 
more about renormalized solutions we refer to [2] and [T5] . 

The usual proof of the renormalization property relies on a commutator 
estimate. It is this estimate that we have to quantify somehow here. More 
precisely we try to bound quantities like 

Wifc= / n IlX ~f\, d \n e (t,x)-n e (t,y)rdxdy, (1.8) 
P ' Jm.™ (\x -y\ + h) d 

uniformly in h. Those norms can be seen as a generalization of usual Sobolev 
norm, in particular we recall that 

I|x-y|<l 



R2d \x-y 



d+2s 



\n e (t,x) - n e (t,y)\ dxdy 



is equivalent to the usual H s norm for s G]0, 1[. This is wrong though for 
s = 0, i.e. \\. 1 1 2,0 is actually stronger than L 2 . In this case p = 2, it is in fact 
easy to see in Fourier that ||.||2,o more or less controls the log of a derivative 
and thus provides compactness. 

We can prove explicit estimates for the norms ( 11. 81) 

Theorem 1.2 Assume (II. 5ft . ( 11.61) . ( 11.71) . that is uniformly bounded in 
L 1 PI L oc (M d ) and is compact in L 1 (M rf ). 3C > only depending on the 
uniform bounds in e s.t. the solution n E (t,x) to ( II. 4p satisfies for any t <T 

' \n £ (t,x) - n £ (t,y)\dxdy 



(\x -y\ + h) d 
I 




\x-y\<2 



(\x -y\ + h) d 



\n° £ (x) — n° £ (y) \ dxdy 



(\x — y\ + h) c 

+ C EL + C \i 0gh \vp 

where p = min(2,p) and 1/p* + 1/p = 



'' '' ' \d e (s,x) — d £ (s,y)\ dxdyds 
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Remarks. 

1. Lemma 13.11 below shows that Theorem 11.21 in fact implies Theorem 11.11 
but its proof is of course more complicated. 

2. In addition of providing an explicit rate, Theorem 11.21 does not require 
the compactness of the sequence a £ . Of course as it is uniformly in L^°W^ ,P , 
it is always compact in space but not necessarily in time. 

3. It is possible to replace (jl.5p by 



The estimate then uses the exponential of this quantity instead of e . 

4. If the sequence Va e is equiintegrable then some kind of rate can also be 
obtained. 

5. Assumption (jl.7p can also be extended by asking r £ to satisfy only 



The norms defined by (11.81) are in fact critical for the problem (II. ID . 
Indeed (II .ip contains the case of the linear transport equation (take / = Id). 
In this last case, one may use the characteristics and it was proved in [IT] that 
one indeed propagates a sort of log of derivative on them. If n° e W 1,p then 
this implies a result like Th. 11.21 Moreover at the level of the characteristics, 
it is not complicated to obtain examples showing that this logarithmic gain 
is the best one can hope for. 

Note that contrary to [11] , we work here at the level of the PDE; because 
of the shocks, the characteristics cannot be used when / is non linear. This 
unfortunately makes the corresponding proof considerably more complicated 
and in particular it forces us to carefully track every cancellation in the 
commutator estimate; we also refer to [5] for an example in a different linear 
situation where a problem of similar nature is found. 

Th. 11.21 gives a rate in | logh] 1 ^ which is probably not optimal. In the 
linear case / = Id, [J2] shows that the optimal rate is 1. In our non linear 
situation, it seems reasonable to conjecture that it should be the same (at 
least for p > 2) but it is obviously a difficult question. 

The proof of Th. 11.21 requires the use of multilinear singular integrals. 
This has been an important field of study in itself (we quote only some results 
below) but quite a few open questions remain, making the optimality of Th. 
11.21 unclear. 




r e \\h,i < C \\n £ \\ hA . 
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The first contributions for multilinear singular integrals were essentially in 
dimension 1, see [7J, [H] or JTU]. The theory was later developed for instance 
in [20], [22], [23]. In dimension 1, an almost complete answer was finally 
given in [27]. In higher dimension, the most complete result that we know 
of, [29], unfortunately does not contain the case that we have to deal with 
here. 

Let us conclude this introduction by mentioning two important and still 
open problems. Of course many technical issues are still unresolved: The 
optimal rate, the case where a £ G BV instead of at least W 1 ' 1 ... 

First of all, in many situations a bound on the divergence of a £ is not 
available. However when / is a logistic function for example, Eq. (II. ip still 
controls the maximal compression, contrary to a linear transport equation. 
It means that this case should actually be easier to handle in the non linear 
setting. 

Second some models do not provide any additional derivative on the ve- 
locity field a. For instance in porous media, one finds the classical coupling 



but one could also consider the non viscous equivalent of (11. 3p . Of course 
the method presented here fails in those cases... 

The next section gives a quick proof of Corollary 11.11 The next section is 
devoted to Th. 11.11 and the last one to Th. 11.21 

In the rest of the paper, C will denote a generic constant, which may 
depend on the time interval [0, T] considered, uniform bounds on the initial 
data n° £ or on a £ but which never depends on e or the parameter h that we 
will introduce. 



Define a sequence of approximations n E , a £ where n £ solves f ll.4p with initial 
data n° and a £ is obtained through n £ by solving ( II. 3p . 
As ( 11. 4p is conservative then one obviously has 



a = 



V0, div x (a(n) V0) = g, 



2 Proof of Corollary 11.1 



n e(t, .)||Li(Rd) 




By the maximum principle 




||w e (t,.)IU 00 CR a ) - ll(/(^) diva e)-IU°°' 
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where (.)_ denotes the negative part. Using (11.21) implies that 

^\\ n e{t, .)IIl°°(r*) < \\{fg{n £ )) „\\ L oo <C||7i e (t,.)|Uo 

by the assumption in Corollary 11.11 Hence by Gronwall's lemma, the se- 
quence n £ is uniformly bounded in L°°([0, T], L 1 PI L°°(M> d )) for any T > 0. 

Thanks to g(0) = 0, the usual estimate for (ll.2p then gives that a £ is 
uniformly in L°°([0, T], W 1 ' p (M d )) for any 1 < p < oo. (Oft . fTL6|) . fTTTjl are 
hence obviously satisfied. 

To apply Th. 11.11 it only remains to obtain the compactness of a £ (note 
that the refined Th. 11.21 does not require it). First we need an additional 
bound on n e . Multiplying Eq. (II. 4p by n e and integrating, one finds 

,0/„,\|2 j„ , 1 / / i |2 



£ / / |Vn £ | dx < I \n (x)\ dx H — / / |n e | diva £ (ia;. 

Thus the previous bounds show that e 1//2 Vn e is uniformly bounded in L 2 . 
Now using the transport equation ( j 1.4ft and the relation (II. 2ft implies for 

*' = /Y 

<9*a £ =V A' 1 d t (g{n £ (t, x))) = -VA^div (a £ /i(n £ )) 

- VA-^g'f - h)(n £ ) diva £ + eVg(n £ ) - eVA' 1 g'\n £ ) \Vn £ \ 2 . 

This proves that dta £ is uniformly bounded in L 2 ([0, T]) with values in some 
negative Sobolev space. Therefore a £ is locally compact in L p (IR d ) with p 
large enough, more precisely p > (1 — l/d)^ 1 by Sobolev embeddings. 

It only remains to control the behaviour at oo of n £ and hence a £ . By 
De la Vallee Poussin, since n° G L 1 , there exists ip £ convex with 

ip(x) — > +oo as |x| —> +oo, V 1 ^ G and s.t. 

ijj(x) \n°(x)\ dx < oo. 

By the convexity of ip, one obtains 

/ |n e (t, ar)| da; < / |VV>| |diva £ | |/(ra e )| dx < C. 

dt J^d J^d 

This implies that Vi G [0, T] 

^(x) |n £ (t,x)|dx < C. (2.1) 
8 



By (11.21) . it has for first consequence that a £ is globally compact in L p , (1 — 
l/d)^ 1 < p < oo. Applying Th. one deduces that n £ is locally compact 
in L l and by (12. ip that n £ is compact in L 1 and so in any L p , 1 < p < oo. 
Let us now extract two converging subsequences (still denoted by e) 



We may now easily pass to the limit in every term of ( II .4p and (II. 2p to deduce 
that n and a are solutions, in the sense of distributions, to ( II .ip coupled with 



Proving that n is an entropy solution to (II. ip follows the usual procedure. 
For any <fi G C 2 convex, we first note that 

d t (j){n e ) + div x (a £ q(n £ )) + (<p'(n £ )f(n £ ) - q(n £ )) div^ a £ < 0, 

with q' = </>'/'. With the compactness of n £ , one may pass to the limit in 
each term and obtain the same property for n, which concludes the proof of 
Corollary 11.11 



3.1 The compactness criterion 

We first introduce the compactness criterion that we use. Define a family 
K h {x) = l/OI + h) d for |x| < 1 and K h non negative, independent of h, 
with support in 5(0, 2) and in C QO (R d \ B(0, 1)). 

Lemma 3.1 A sequence of functions Uk, uniformly bounded in L p (M. d ) is 
compact in Lf oc if 

limsup | log | 1 / K h (x — y) \uk(x) — Uk(y)\ p dx dy — > as h — > 0. 

k jR2d 

Conversely if Uk is globally compact in L p then the previous limit holds. 
Proof. We recall that Uk is compact in LP iff 




3 Proof of Theorem 11.11 




as rj — > 0. 
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So assuming u k is compact, one simply decomposes 

sup/ K h (x - y)\u k (x) - u k (y)\ p dxdy < C 



k 



+ C V sup/ 2 d > fc (z)- Wjfc ( 2 /)p>d ; zeft/ 

n<|iogh| k h ' n ~ l >\*-y\< 2 ~ n 

<c + c <K 2 ~ n )> 

n<\ log h 

which gives the result. 
Conversely assume that 

a(h) = limsup | log h\^ x I Kh(x-y) \uk(x) — u k (y)\ p dx dy — y as h — > 0. 

k jR2d 

Denote Kh{x) = Ch\ \ogh\~ 1 Kh(x — y), with Ch s.t. 

Kh(x) dx = 1, 



and therefore Kh a convolution kernel. Note that Ch is bounded from below 
and from above uniformly in h. Now 

\\u k - K h * x u k \\ p LP <\ \ogh\~ p / (/ ^(ar-^M^-u^ldy <&/ 

<| \ogh\~ p WKhW' 1 ^ 1 / K h (x - y)\u k (x) - u k (y)\ p dydx 

JK 2d 



<C\\ogh\ / K h (x - y)\u k (x) - u k (y)\ p dy dx 

is converging to uniformly in k as the limsup is and it is converging for 
any fixed k by the usual approximation by convolution in L p . On the other 
hand for a fixed h, Kh * u k is compact in k and this proves that u k also is. 

3.2 The main argument given for a linear transport 
equation 



Before proving Theorem II. 1[ we wish to explain the main idea behind the 
proof in a simple and wellknown setting. Let us consider a sequence u £ of 
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solutions to the transport equation 

d t u £ (t, x)+v £ - Vu £ (t, x) = 0, t e [0, T], x e R d , 

U £ (t = 0, x) = Ue(x), 

for a given velocity field. The following result was originally proved in [T8] 

Theorem 3.1 Assume that vP e is uniformly bounded in L 1 flL 00 and compact. 
Assume moreover that v e is compact in LP ', uniformly bounded in L^°W^' P for 
some p > 1 and that divt> £ = 0. Then the sequence of solutions u £ to (13.1 j) 
is compact in L 1 . 

Proof of Theorem 13.11 

First of all notice that u £ is uniformly bounded in L ( °°(L],nL™). Moreover 
as v £ is compact, one may freely assume that it converges toward a limit 
v E L^°W^' P (by extracting a subsequence). 

Now define 

Q £ (t) = / K h (x - y) \u £ (t,x) - u £ (t,y)\ 2 dxdy. 

From Equation (13. ip the divergence free condition on v £ , one simply computes 

— - = / VK h (x -y)(v e (t,y) -v e (t,x))\u e {t,x) - u £ {t,y)\ 2 dx dy. 
at j R 2d 

Therefore by introducing the limit v 

dQ 



<C \\v £ - v\\ L p \\VK h \\ L i 



dt 

+ / VK h (x-y)(v{t,y)-v(t,x))\u e {t,x) - u £ {t,y)\' 2 dx dy. 



The second term is equal to 



/ / (x-y)®VK h (x-y) : Vv(t,9x+ (1 -8)y) \u £ (t,x) -u £ (t,y)\ 2 dxdy, 

JO JR 2d 

with A : B denoting the full contraction of the two matrices. Note that for 
\x\ > 1, VKh is bounded and for \x\ < 1, 



___ f \ XtSX 

x g> VK h (x) 



{\x\ + h) d+1 \x[ 
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Define 

~ _ Id \x\ 

K h (x) = x ® VK h {x) - A (N+ '^ d+1 I M < ls 

with A = f S d-i w\ du. 

Thanks to the definition of A, Kh is now a Calderon-Zygmund operator, 
meaning that for any 1 < q < oo, there exists a constant C independent of 
h s.t. 

||^*fflU«<cy| L ,. 

As v is divergence free, one may simply replace by K h 

[ [ (x-y)®VK h {x-y):Vv{t,9x + (l-9)y)\u £ {t,x)-u £ {t,y)\ 2 dxdy 
Jo Jm 2d 

= I I K h (x-y):Vv(t,9x + (l-9)y)\u £ (t,x)-u £ (t,y)\ 2 dxdy 
Jo Jm. 2d 



< C [ [ K h (x-y) \Vv(t,0x+ (1 -9)y) - Vv(t,x)\dxdy 

Jo JR 2d 

+ \\K h *(Vvul)\\ L i+2\\u e K h *(Vvu e )\\ L i + \\u 2 e k h *Vv\\ L i. 

Thanks to the uniform bounds on u £ , and changing variables, one immedi- 
ately deduce that 

/ / {x-y)®VK h (x-y): Vv(t,9x + (l-9)y) \u £ (t,x) - u £ (t,y)\ 2 dx dy 
Jo Jm 2d 

<C + C [ K h {x-y)\Vv{t,x)-Vv{t,y)\dxdy. 

jR 2d 

Putting together all the terms in the estimate, we have 

dt h 
+ C K h (x-y)\Vv(t,x)-Vv(t,y)\dxdy 

jR 2d 



or 

'\v e - v\\ L p 



Qe(t) <C + C- 

+ C f f K h {x-y)\Vv{t,x)-Vv{t,y)\dxdy 

Jo JR 2d 



h 

t-T 



+ C f K h (x-y)\n £ (x)-n° £ (y)\ 2 dxdy. 

jR 2d 
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As n° £ is compact and v is independent of e then the previous estimate shows 
that 

lim | \ogh\~ l limsupsup / K^(x — y) \n £ (t, x) — n £ (t, y)\ 2 dx dy = 0. 
h ->° e t Jm 2d 

Lemma [3.11 then proves that u £ is compact in space. However by Eq. ( 13. ip . 
dtu £ is uniformly bounded in L™(W~ l,p ). Therefore compactness in time 
follows and the theorem is proved. 

3.3 A simple proof for Theorem 11.11 

We first give here a simple proof of the compactness. This proof is not 
optimal in the sense that it does not give an explicit rate for how the norm 
in our compactness criterion behaves 

T f 

K h {x - y)\n £ (t,x) - n £ (t,y)\dxdydt. 



o 

This is however a more difficult problem, which is partially dealt with in the 
next section. 

As a £ is compact in L p , by extracting a subsequence (still denoted by e), 
a £ converges strongly in L p to some a G W 1,p . By the compactness of d £ 
and n° and by Lemma 13.11 we may assume without loss of generality that 
there exists a continuous function 8{h) with 5(0) = 0, independent of e and 
a function ot(e), s.t. 



T f 

/ K h (x — y) \d £ (t, y) — d £ (t, x) \ dx dy dt < 8(h), 
°/f (3.2) 



\a £ (t,x) — a(t, x)\ p dxdt < a p (e). 

Note that the estimate is written for Va and not for the sequence Va £ as no 
compactness can be assumed on Va e . 
Then one proves 



log 


h\- 


log 


h\- 


log 


h\- 


f 




h 


s 
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Proposition 3.1 Let n £ be a sequence of solutions to (II. 4p with initial data 
n° £ uniformly bounded in L 1 fl L°° and compact in L 1 . Assume (ll.5p . (II. 6p . 
(11. 7p and hence (13. 2p . Then for some constant C uniform in h and e 



The disappointing part of Prop. 13.11 is that the rates 5(h) and a(e) are not 
explicit but depend intrinsically on the sequence a e . See the next section for 
a more explicit (but much more complicated) result. 

Prop. 13. II proves the compactness in space of n £ by Lemma I3TT1 The com- 
pactness in time is then straightforward since n £ solves a transport equation 



Hence Theorem 11.11 follows. 
Proof of Prop. HO . 

The proof mostly follows the steps of the proof of Theorem 13. II The main 
differences are the nonlinear flux, the vanishing viscosity terms and the fact 
that now the field a e is not assumed to be divergence free (only bounded). 

First of all, by condition (II. 6p . for any T > 0, n e (t,x) is bounded in 
L 1 n L°°([0, T] x R d ), uniformly in e. 

We start with Kruzkov's usual argument of doubling of variable. If n £ is 
a solution to (II .4p then 

d t \n £ (t,x) - n £ (t,y)\ + div x (a £ (t, x)F(n £ (t, x), n e (t, y))) 

+ div y (a £ (t,y)F(n £ (t,y),n £ (t,x))) + div x a £ (t,x) G(n £ (t,x),n £ (t,y)) 

+ div y a £ (t,y) G(n £ (t,y),n £ (t,x)) - e 2 (A x + A y ) \n £ (t,x) - n £ (t,y)\ < 0. 

This computation is formal but can easily be made rigourous by using a 
suitable regularisation of |.|. Here F satisfies 




K h (x—y) \n £ (t,x)—n £ (t,y)\ dxdydt < C -r^+C 5(h) \ \ogh\+C 




F'(C,0 



/'(Osign(f-C), F(Z,() = 0, 



which means that 



F(C = (/(0 - /(C)) sign(e - C) = F((, 0- 



And as for G 



G((,0- 
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Now define 

Q(t) = / K h (x - y) \n £ (t,x) - n £ (t,y)\ dxdy. 

Remark that 

K h (x - y) (A x + A y ) \n £ (t,x) - n £ (t,y)\ dxdy 

i 

= e 2 AK h (x — y) \n £ (t,x) — n £ (t,y) \ dxdy 



e 1 



<Ce 2 \\AK h \\ L ,<C^-. 
Using this and because of the symmetry of F and the antisymmetry of G 

Jt Q{t) ~ C h + (Ix-yf+V 2 ' ^ V) ~ 

F(n £ (t,y),n £ (t,x)) dxdy 

+ / VK(x - y) ■ (a £ (t, y) - a e (t, x)) F(n £ (t, y),n £ (t, x)) dx dy 

J\x-y\>l 

+ / K h (x - y) (diva £ (t,x) - div a e (t,y)) G(n E (t,y),n e (t,x)) dxdy 



e 2 

--C— + A + B + D. 



Let us begin with the last term. Use ( 11. 7ft to decompose 

D < / K h (x - y) \d £ (t,x) - d e (t,y)\ \G{n s (t,y),n e (t,x))\ dxdy 



+ C K h (x-y)\n £ (t,x) -n e (t,y)\ \G(n E (t,y),n e (t,x))\dx dy. 
Jm. 2d 

As G(n E (t,x),n £ (t,y)) is uniformly bounded in L°°, one gets from (13. 2p 

I Ddt<\\ogh\5{h) + C [ Q{t)dt. (3.3) 
Jo Jo 

For the second term B, just note that VK 6 C™(R d \ B(0, 1)), and that 
\F(n e (t,x),n £ (t,x))\ < \f(n £ (t,x))\ + \ f(n £ (t,y))\ is uniformly bounded in 
L l . So one simply has 

B < C. 
15 



The main term is hence A. Using again (13. 2p and the bound on \F(., .)|, one 
gets 



a(e) / / x — y 



Ait< - c ^ + L 7^ ffl <i*-,i + *)*M*-»r (a(t - B) - a( ^ 

F(n £ (t,y),n £ (t,x)) dxdy dt 
<_C a M + ff{ ( '7»f< J - g ) , =Va(t, fa + (!-%) 

ft io 7o y|*-»i<i (f - y| + ft) + k - y\ 

F(n £ (t, y),n £ (t, x)) dx dy d6 dt. 

Still using fl3~2|) . 

T Adt<C f^ + |logft|5(ft) 

+ / / n} X _~. V ^}Z~\^ ~ ■■Va(t,x)F(n e (t,y),n e (t,x))dxdydt 
a(e 



o J\x-y\<i {\x - y\ + h) d+1 \x - y\ 

T 

I Ml i \ I 

< c 



+ \\ogh\8{h) S j + J E{t)dt. 









(1*1 





Denote as in the proof of Theorem 13.11 

A / u(dS(u), K h (x) 1 

By the definition of A, is a Calderon-Zygmund operator bounded on any 
LP for 1 < p < oo. Now write 

F — K h (x - y)Va(x) F(n £ (t,y),n £ (t,x)) dxdy 

JR 2d 

r | ^ — ^ | 
+ A / j —— z —diva(t,x)F(n £ (t,y),n £ (t,x))dxdy 

J\x-y\<l \\ x ~ V\ + n ) 

< / K h (x-y)Va(x)F(n £ (t,y),n £ (t,x))dxdy + CQ(t), 



as the divergence of a is bounded. 
Introduce 

Xe(t,X,Q = h<Z<n s (t,x)- 
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Then note that Xe is compactly supported in £ and that 



F(n £ (t,y),n £ (t,x)) 



POD 

/ f'(0\Xe{t,X,0-Xe{t,y,0\ 2 ^. 

Jo 



Hence as Va G L p , and Xe is uniformly bounded in Lf^(L\ fl for l/p + 
l/p* = l, 

Kh(x — y) Va(x) F(n £ (t, y),n £ (t, x)) dx dy 

i 

= f'(0 / K h {x - y) Va(x) |x e (i,x,0 - Xe(t,y,0\ 2 d £ dxd V 

< [ / , (0(ll^fc*(Vox^)|Ui + 11^*^)11^ +2||^ fc *(Vox e )||L 1 )^ 

< c. 

Combining all estimates we conclude that 

Q(t) <Q(0) + C^ + C\\ogh\6(h) + C^ + C [ Q(s)ds. 

h h Jo 

The initial data Q(0) is bounded by ( 13. 2h and finally by Gronwall lemma we 
obtain on any finite interval 



e 2 „ a(e) 



Q{t) < C — + C | \ogh\ 5(h) + C ■ 



h 2 1 & 1 w /i ' 

which proves the proposition. 

4 An explicit estimate : Proof of Theorem 



Checking carefully the proof of Prop. 13.11 one sees that to get an explicit 
rate, it would be necessary to bound a term like 

VK h (x - y) (a £ (x) - a £ (y)) \g £ (x) - g £ (y)\ 2 dxdy (4.1) 

only in terms of the W 1,p norm of a £ and the L l fl L°° norms of g £ . 

Here we do not aim at optimal estimates, just explicit ones. We present 
a very elementary proof of 
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Proposition 4.1 Let 1 < p < oo, 3C P < oo s.t. Va(x), g>(:c) smooth and 
compactly supported 

VK h (x - y) (a(x) - a(y)) \g(x) - g(y)\ 2 dxdy 

i 

< C p \\g\\ L ™ \\g\\&niy l|Va|| L p nL i \ \ogh\ 1/p 

+ C p (||diva|| L oo + ||Va|| LP ) / K h (x - y)\g{x) - g{y)\ 2 dxdy, 

with 1/p* + 1/p — 1 and p = min(p, 2) . 

Note that the rate | logh] 1 ^ is most probably not optimal. A way to obtain 
a better rate could be to combine Lemma 14.11 below with the estimates in 
[29] as we suggest below. 

The kind of Calderon-like estimate like Prop. 14.11 has been extensively 
studied in dimension 1, see for instance [7] or [9], [10]. The situation in higher 
dimension is however more complicated. In particular it seems necessary to 
use the bound on the divergence of a e to estimate (14.1 1) (as was already 
suggested by the proof of Prop. 13. ip . 

Following the previous section, a simple idea would be to estimate (14. ip 

by 

C ||diva|| L oo / K h (x -y)\g £ (x) - g £ (y)\ 2 dxdy 



+ [ [ L h {x-y) :Va{6x + {l-6)y)\g e {x) - g £ (y)\ 2 dx dy, 
Jo Jm. 2d 

where Lh is now a Calderon-Zygmund operator. Expanding the square, one 
sees that it would be enough to bound in some L q space 

/ [ L h (x-y):Va{dx + {l-d)y)g £ (y)dy. 

JO JR d 

Using Fourier transform (we denote by T the Fourier transform) and an easy 
change of variable, this term is equal to 

jR 2d 

with 
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We now have a multi-linear operator in dimension d of the kind studied in 
Muscalu, Tao, Thiele [2S]. Unfortunately m does not satisfy the assumptions 
of this last article as it does not have the right behaviour on the subspace 
£1 || £2- Instead it would be necessary to have a multi-dimensional equivalent 
of [27] (which, as far as we know, is not yet proved) or to use Lemma 14.11 

4.1 Proof of Theorem 11.21 given Prop. 14.11 

For the moment let us assume Prop. 14.11 Define 

Q(t) = / K h (x - y) \n e (t,x) - n e (t,y)\ dxdy. 

We follow the same first steps as in the proof of Prop. 13. 1\ with the same 
notations. We obtain 

<1Q <C + C^ + CQ(t)+C [ K h (x-y)\d e (t,x)-d e (t,y)\dxdy 



dt ~ h 2 

+ / VK h (x -y) ■ {a £ {t,x) - a e {t,y)) F(n e (t,y),n e {t,x))dxdy. 

jR 2d 

We only have to bound the last term. Let us introduce again 

Xe(t,X,£) = h<£<n,(t,x)- 

Note that Xe is supported in £ in [0, H^Hl 00 ] C [0, C]. 
Now write 

VK h {x - y) ■ (a e (t, x) - a £ (t, y)) F(n E (t, y), n E (t, x)) dx dy 

/'(C) / VK h {x - y) ■ {a e {t, x) - a e {t, y)) 

\Xs(t,y,Q - Xs{t,x,€)\ 2 dxdyd£ 

rC 







<C7|log/i| 1 - 2 ^ + C / K h (x-y) [ \xs(t,y,0-Xs(t,x,0\ 2 ^dxdy, 

Jm. 2d JO 



using Prop. 14.11 and the uniform bounds on ||a E ||L°°iS anc ^ IIXelUinL 00 - Now 
simply note that because of the definition of Xe 



c 

2 







Xe(t,y,0 - Xe(t,x,£)\ < \n £ (t,x) - n E (t,y)\, 
19 



and the last term in the previous inequality is hence simply bounded by Q. 
One finally obtains 

^r<C + +C^ + C Q(t) + C | log hi 1 - 2 /? 
at h z 

+ C K h (x -y)\d £ (t,x) - d E (t,y)\dxdy. 



To conclude the proof of Theorem II .2\ it is now enough to apply Gronwall's 
lemma. 



4.2 Beginning of the proof of Prop. 14.11 

As before we will control a(x) — a(y) with Va. Contrary to the previous 
case though, it is not enough to integrate over the segment. Instead use the 
lemma 



Lemma 4.1 

f ( x — y \ dz 
a,i(x) - cii(y) = \x - y\ I ip[z,- ■ Va { (x + \x- y\ z) 

</b(o,i) V \ x ~y\J \ z \ 

f . ( x — y \ ^ , . . dz 

+ \x-y\ / ip[z,- , ■ Va, (y + \x - y\ z) — -— r , 

Jb(o,i) V \x-y\J \A 

where \z\ip is Lipschitz on B(0, 1) x S and for a given constant a, 

( x-y\ dz x-y 



b(o,i) v \x-y\j \ z \ \ x ~y 



Proof of Lemma 14.11 We refer to [8] for a complete, detailed proof. Let us 
simply mention that the idea is to integrate along many trajectories between 
x and y instead of just the segment. □ 

Lemma 14.11 gives two terms that are completely symmetric and it is 
enough to deal with one of them. After an easy change of variable, one 
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finds 

VK h (x - y) (a(x) - a(y)) \g(x) - g(y)\ 2 dxdy 



/ / 7 r , \ j i -i / / ip(z,u)®u : Vafx + rz) 

dz 

\g(x) — + ru) \ 2 dtu _ 1 dr d:r + symmetric, 



where A : B denotes the total contraction of two matrices Yliij ^ij Bij. 
Now define 

L(z, oj) = ip(z, u>) <S> ou — A Id, 

Note that 

/ / 7 r , yn / / : Va(a; + rz) 

y Rd y {r + h) d + i j mi) j sd -r y 



dz 

\g{x) — g(x + ru)\ 2 dbj dr dx 



Z 1 /" r d L(z,lo) _ . . . , , , mo 1 

< / / : va(i + rzj |<7(x) — + rwj 1 



/o iRdxB(o,i)xSd-i (r + a) d+1 ' " '' ' ' ' \z\ d 1 

+ C||diva|| L oc f [ ^ \ _ 2 

By the definition of .K^, the second term is bounded by 

C ||diva|| L oc / K h (x - y)\g(x) - g(y)\ 2 dxdy. 

and it only remains to bound the first one. In order to get the optimal rate 
for Va G L p with p > 2, we need to introduce an additional decomposition 
of Va. For p > 2 as W may be obtained by interpolating between L 2 and 
let 

Va = A + A, p|| £00 < 2||Va|| LP , \\A\\ L 2 < 2||Va|| LP . 

If p < 2 then we simply put A = Va. In both cases, if Va is smooth and 
compactly supported then one may of course assume the same of A and A. 
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Define 

r d L(z,u) Al s t s j dz 

'0 JB(0,l)xS d ~ 



Q(A, g) = I I ^ - : A(x + rz) g(x + ru) du dr. 



The term with A may be bounded directly by using the L°° norm of A; for 
the other one simply by expanding the square \g(x) — g(y)\ 2 , one obtains 

VK h (x-y) (a(x) - a(y)) \g(x) - g(y)\ 2 dxdy 

i 

< C (\\diva\\ L °° + ||Va|| L p) / K h (x - y) \g(x) - g(y)\ 2 dx dy 
+ [ (-2gQ(A,g)+g 2 Q(A,l))dx. 

Note that bounding Q(A, 1) is in fact easy as it is an ordinary convolution and 
- L defines a Calderon-Zygmund operator. However the control of Q(A,g) 
essentially requires to rework Calderon-Zygmund theory. 
Of course for r of order h then one has 

* f r d L(z,u) dz 

/ „ (r + hV+i : A (x + rz)g(x + ru)du-- rT dr 

JB(Q,l)xS d - 1 V + n ) \ z \ L i 

< — f f \\A(x + rz) g(x + ru)\\ L1 du dr 
h JO J BiO^xS 11 - 1 " \ z \ 

< \\A\\ L i \\g\\L°°- 

It is hence enough to consider 

p\i a \ f 1 f r d L(z,u) ., . , s , dz , 

Q(A, g)= / r-T^i = A(x + rz) g(x + ru) du j-^ dr. 

Jh JB(0,l)xS' d - 1 V + n ) \ z \ 

We introduce the Littlewood-Paley decomposition of A (see for instance 
Triebel [31]) 



A(x) = ^ Ms), 



i=0 



where for i > 0, A\ = J r Ap(2~ 1 ^) with p compactly supported in the annulus 
of radii 1/2, 2; and A = TAp (£). The functions p and p determines 
a partition of unity. We note either Tg or g the Fourier transform of a 



22 



function g. In the following we denote by Vi the projection operator Vi4> = 
T~ x v{ 2^00- 

There is an obvious critical scale in the decomposition which is where 1~ % 
is of order r. Accordingly we decompose further 

Q(A,g) = Q 1 (A,g) + Q 2 (A,g) 

rt- 1 r 



V- f f r d L(z,u) . . . . . J dz , 

Ef l f r d L(z,u) . , . . , r dz r 

/ / t — — 7— r : + rz) olx + ran aa> — ; — -dr. 



Each term is bounded in a different way. Note of course that in Qi as r > h 
there is of course no frequency % higher than | log/i| (they are all in Q 2 ). 

4.3 Control on Q\ in L 2 

The aim is here is to prove 

Lemma 4.2 VI < q < oo, 3C > such that for any A and g smooth and 
compactly supported functions, 

||Q 1 (As)IU 1 <C'||A|| s o i \\g\\ Lq *. 

where B® x is the usual Besov space and 1/q* + 1/q = 1. 

As we wish to remain as elementary as possible here, we avoid the use of 
Besov spaces in the sequel. Instead for q = 2 it it possible to obtain directly 
the Lebesgue space by losing | log/i] 1 / 2 namely 

Lemma 4.3 3C > such that for any A and g smooth and compactly sup- 
ported functions, 

\\QM,g)\\ L i < C | log/il 1 / 2 NsrlU- \\A\\ L2 . 



The proof is relatively simple. Indeed in Qi since r < 2 l , Ai does not 
change much over a ball of radius r. Therefore, we simply replace A^x + rz) 
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by Ai(x) in Qi- This gives 
Qi(i4,0)</ + // 

i<\\o g h\ Jh Jb<p,i)xS*-i {r + h) \z\ 

ft' 1 i r (4.2) 
+ E / ~Th / |^ + w)|dw 



i<\logh\ 



/ \Ai(x + rz) - Ai(x) 
Jb(o,i) 



dz 

■ dr. 



\z 



d-l 



Let us bound the first term. As Ai does not depend on z anymore, this term 
is simply equal to 



Ai(x) / / Li(ruj)gj(x + ru) r d 1 dudr, 

i,j<i 



where 



t I \ r I \ t\ f T 'K<2- i T I \ ^ 

By the definition of A, A = j sd -i w\ doo and hence Lj is a Calderon-Zygmund 
operator with operator norm bounded uniformly in %. 
Now write for 1/q* + 1/q — 1 

l|J(*)IUi< E IIA||L.||^*ff||^ 

i<| log h\ 

C\\g\W Y.W A i\\L« = C\\g\\ Lq * \\A\\ B o y 



< 



Let us turn to the second term. 



\\H\W< E / / '-^r-l \M* + rz)-Mx 
i<^ gh \ Jsd - lJh r + h ^(0,1) 

<C\\g\\ Lq * f f \\M- + rz)-A 



d-l 



L 1 
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So 



\n\W<c\\ g \\ Lq * (IIA(.)IU. + I|A + i(.)IU ? + IIA-i(-)IIlO 

i<| log /i 



2 r Lz 



r + /i \z 



d-i 



'B(0,1) ./h 

where we used the localization in Fourier space of the Aj and more precisely 
the well known property 

\\M- + v)- M-)\\l« < ct \n\ (|K(.)IU« + II^WOIU* + 

One then concludes that 

\\II\W<C\\g\\ Lq * WMl* = C\\g\\ Lq * \\A\\ B o y 

i< | log h\ 



Combining the estimates on / and II in (14.21) gives Lemma 14.21 

For the proof of Lemma 14.31 it is enough to observe that in the case q = 2 
by Cauchy-Schwartz 



ii^iu^iiog/^ hrp 



] =|log/ i | 1 / 2 p|| L2 . 



i< | log h\ 



4.4 Control on Q 2 for A G L 2 

As for usual Calderon-Zygmund theory, the optimal bound on Q 2 is obtained 
in a L 2 setting namely 

Lemma 4.4 3C > s.t /or any g and A smooth with compact support 

\\Q2(A,9)\\v<C\\g\\ L » \\A\\ L2 . 



To prove this, first bound 

\Q2(A,g)\ < \\g\\ Loo [ [ 

Jo Js*- 1 



r d Liz, uj) 



B(0 



,i) + h) 



d+l ' 



dz 



Mx + rz)--^ 

i>-log 2 r 



duj dr. 
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Hence 



\\Q2{A,g)\\»<CR\\g\\l a + £ f f 

K Jo Js*- 1 



r d L(z,u) 
(r + h) d + l ' 



B(0,1) 



dz 



rz) 



Id— l 



i>— logo T 



du dr. 



L-' 



Use Fourier transform and Plancherel equality on the last term to bound it 
by 



/ r 4- 



a,P= 



1 Jo r + h Js d ~ 1 J\z\>i/ 



'B(0,l)xB(0,l) 

One only has to bound the multiplier 



L a p(z,u)L a p(z',u))e i t- r ( z z>) 7-^=1 i Md-i d £ dudr 



\zr ^ \z' 



m(£ } u,r) 
Define 



LcftMLaeV^e**'-* ^ <hJ 



fi(0,l)xB(0,l) 



Z I ^ ^ 1 2;' I d—1 ' 



M(£,u,r,s) = / L a p(su,uj)e irs t- U du, 



S d-i 



such that 



1 r i 




M (£, oo, r, s) M(£, u, r, s') ds ds' . 



JO 



Assuming for instance that £ is along the first axis, by the regularity on ip 
and hence L given by Lemma 14.11 



\M(£,u,r, s) 



< 



S d-i 
C 



L aP (su,u)e irs ^ Ul du 



< 



r\Z\ 



S d-i 



\V z L a/3 (su,oj)\ du 



rs\£\' 



As M is also obviously bounded, one deduces that 

C 



\M(£,u,r,s)\ < 
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Introducing this in m immediately gives 

C 



r If I 

Therefore eventually 



\\Q2(Ag)\\ L2 <CR\\g\\l ao + ^ f l^(OI 2 f 

K £fL x j^ d J\t\- x r + ft r Is I 

3 

5^ / 1^(0 1 2 « 

a,fl=l- /Rd 



«)/8 



< c IMU°° II^IU 2 ; 

by optimizing in R, which proves the lemma. 
4.5 Control on Q for A G L p 

To get an optimal bound, one should now try to obtain weak-type estimates 
on Q 2 , showing for instance that it belongs to L 1 — weak if A e L 1 ; and then 
use interpolation. Additionally, we would have to use the bound given by 
Lemma 14.21 with Besov spaces. 

However here we will be satisfied with any explicit rate, even if it is 
not optimal. We hence completely avoid some (not negligible) technical 
difficulties and obtain instead 

Lemma 4.5 VI < q < oo, 3C > s.t. for any smooth g and A with compact 
support 

\\Q{A,g)\\ L i +L9 < C \logh\^ \\g\\ L oo nL 2 \\A\\ Lg , 
where q = min(g,g*) with 1/q* + l/q = 1. 

Remark. Note that thanks to our decomposition of Va, we only use Lemma 
1431 for q < 2. 

Proof of Lemma 14.51 Fix g and consider Q(A, g) as a linear operator on 
A. The easy control for r < h, Lemmas 14.31 and 14.41 imply that this operator 
is bounded from L 2 to L 1 + L 2 with norm C (\\g\\L°° + (log/i) 1 / 2 ||g|| L 2). 
On the other hand, one has the easy estimate 

* x 1 f dz 



\Q(A,g)\<C\\g\\ L - [ —J— / \A(x + rz 
Jo r + hj mi) 



d-1 
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Therefore for any 1 < q < oo, Q is bounded on L q with norm less than 
C \\g\\L°° I log/i|. By usual interpolation, one deduces the lemma. 

4.6 Conclusion on the proof of Prop. 14.11 

By subsection 14.21 



< C (||diva||ioc + ||Va||xp) / K h (x - y)\g(x) - g(y)\ 2 dxdy 

Jm. 2d 

+ 2||/|| LP . ||g(Ai)ll LP- 

Bound directly Q(A,g) by Lemma [4.51 and observe that Q(A, 1) is bounded 
on any LP with 1 < p < 00. This completes the proof of the proposition. 

Acknowledgement. P-E Jabin is grateful to P.A. Markowich for having 
pointed out this interesting problem. 
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